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O , Abstract. Progress on nuclear liquid gas phase transition (LGPT) or critical 

' behavior has been simply reviewed and some signals of LGPT in heavy ion collisions, 

i especially in NIMROD data, are focused. These signals include the power-law charge 

distribution, the largest fluctuation of the fragment observables, the nuclear Zipf law, 

caloric curve and critical exponent analysis etc. 

o : 

^ ! 1. Introduction 

: 

Phase transition and critical phenomenon are extensively debatable subjects in the 
natural sciences. Recently, the same concept was introduced into the astronomical 
objects P as well as the microscopic systems, such as in atomic cluster [2] and nuclei 
Nucleus is a finite-size droplet-like quantum system in the ground state, it may 
experience the liquid-gas phase transition or critical behavior when the nucleus is 
heated. Till date, various experimental evidences have been accumulated to corroborate 
^ ■ the nuclear liquid gas phase transition or critical behavior. For instance, violent 
heavy-ion collisions break the nuclei into several intermediate mass fragments (IMF) 
instantaneously |H El EI, which might be viewed as a critical phenomenon as 
observed in fluid, atomic and other systems. The sudden opening of the nuclear 
multifragmentation and vaporization [B] channels can be interpreted as the signature 
of the boundaries of phase mixture. In addition, temperature plateau of the nuclear 
caloric curve [9, in a certain excitation energy range gives a possible indication of phase 
transition [TU1 ITT) IT^l [T5] . On the other hand, the extraction of critical exponents in 
the charge or mass distribution of the multifragmentation system ^3] points to the 
continuous phase transition. More recently, the negative microcanonical heat capacity 
was experimentally attempted in nuclear fragmentation which may relate to the 
liquid-gas phase transition (THj, and in atomic cluster [T2j which relates to solid-liquid 
phase transition |18j . respectively. Moreover, the evidence of spinodal decomposition in 
nuclear multifragmentation was recently reported experimentally fHl ED] , which may 
show the presence of liquid-gas phase coexistence region in the finite nuclear systems. 
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Bimodal distributions of some chosen variables measured in nuclear collisions were 
recently proposed as a signature of a first order phase transition in nuclei [21] . A-scaling 
of the largest fragments was also investigated recently and it is probably a probe to 
detect the phase change [221 • The nuclear Zipf's law was proposed to diagnose the onset 
of liquid-gas phase transition or critical behavior [2B1 EH [2Ej and it is experimentally 
supported [231211 • Phase coexistence diagram was also tentatively constructed based 
on the EOS data [2E1I2H1- 

While, extensive theoretical models have been developed to treat such a phase 
transition in the nuclear disassembly, e.g. percolation model, lattice gas model, 
statistical multifragmentation model and molecular dynamics model etc (e.g. see some 
recent articles [3U1 Ell E21 EH1 Ell ES] and references therein). 

In this article which dedicates to Prof Fu-Jia Yang on the occasion of his 70th 
birthday, I will focus on some experimental aspects of the nuclear LGPT or critical 
behavior in heavy ion collisions, specially for NIMROD data. 

2. Signals of the Nuclear Liquid Gas Phase Transition or Critical Behavior 

2.1. Fisher Droplet Model: Power Law of Charge/Mass Distribution 

The Fisher droplet model has been extensively applied to the analysis of 
multifragmentation since the pioneering experiments on high energy proton-nucleus 
collisions by the Purdue group [HIEniEZI- Relative yields of fragments with 3 < Z < 14 
could be well described by a power law dependence A~ T and it was suggested that this 
might reflect the disassembly of a system whose excitation energy was comparable to its 
total binding energy jHEj- The extracted value of power-law exponent was 2 < r < 3, 
which is in a reasonable range for critical behavior 38J. The success of this approach 
suggested that the multi-fragmentation of nuclei might be analogous to a continuous 
phase transition from liquid to gas which is observed in more common fluids. 

In the Fisher Droplet Model the fragment mass distribution can be represented as 



where Yq, r, X and Y are parameters. However, at the critical point X = 1 and Y = 
1, therefore the cluster distribution is given by a pure power law 



The model predicts a critical exponent r ~ 2.21. 

We collaborated an experiment with a 47r-detector combination of charged particle 
ball with neutron ball (NIMROD) in Cyclotron Institute, Texas A&M University to 
explore the quasi-projectile fragmentation from Ar + Al, Ti and Ni at 47 MeV/nucleon 
[2*5] . Well defined QP sources have been reconstructed using a new method which is 
based on three source fits and the Monte Carlo sampling for the assignment of quasi- 
projectile (QP) and their disassembly features have been analyzed. For simplicity, we 
call this experimental data as TAMU data [2E]. The details can be found in Ref. j2E]- 



Y(A) = Y A- T X A2/3 Y A , 




Y(A) = YqA~ t . 



(2) 
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Figure 1. (Left Part) Charge distribution of QP in different E* /A (in unit of 
MeV/nucleon) window. From the TAMU data |2fj) . 

Figure 2. (Right part) The effective Fisher-law parameter (r e //) (a), the effective 
exponential law parameter (A e //) (b), (S2) (c), NVZ fluctuation (d), the mean charge 
number of the second largest fragment (Zi max ) (e), the Zipf-law parameter A (f). See 
details in text. Figure is taken from the TAMU data . 



Fig. presents yield distributions dN/dZ, observed for nine different excitation 
energy (E*/A) intervals for the QP from the reactions of m Ar + 58 Ni from TAMU 
data. At low excitation energy a large Z residue always remains, i.e. the nucleus is 
basically in the liquid phase accompanied by some evaporated light particles. When 
E*/A reaches ~ 5.6 MeV/nucleon, this residue is much less prominent. As E*/A 
continues to increase, the charge distributions become steeper which indicates that the 
system tends to vaporize. To quantitatively pin down the possible phase transition 
point, we use a power law fit to the QP charge distribution in the range of Z = 2 - 7 to 
extract the effective Fisher-law parameter r e jf by dN/dZ ~ Z~ Te ff in Fig. ^ Fig. Efa) 
shows r e ff vs excitation energy, a minimum with r e ff ~ 2.3 is seen in the E* /A range 
of 5 to 6 MeV/nucleon. r e ff ~ 2.3, is close to the critical exponent of the liquid gas 
phase transition universality class as predicted by the Fisher's droplet model 

While, assuming that the heaviest cluster in each event represents the liquid phase, 
we have attempted to isolate the gas phase by event-by-event removal of the heaviest 
cluster from the charge distributions. We found that the resultant distributions are 
better described as exponential form e~ XeffZ . The fitting parameter A e // was derived 
and was plotted against excitation energy in Fig. |2fb). A minimum is seen in the same 
region where r e ff shows a minimum. 
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2.2. The largest fluctuation 

One of the well known characteristics of the systems undergoing a continuous phase 
transition is the occurrence of the largest fluctuations. These large fluctuations in cluster 
size and density of the system arise because of the disappearance of the latent heat at 
the critical point. In macroscopic systems such behavior gives rise to the phenomenon 
of critical opalescence [40J. 

To further explore this region we have investigated other proposed observables 
commonly related to fluctuations and critical behavior. Fig. |2^c) shows the excitation 
function of the mean normalized second moment, (S2), defined as S 2 = M 2 /Mi, where 
M 2 and Mi is the second and first moment. Generally, the k-th moment is defined as 



where m is the multiplicity of the fragment A and the largest fragment A max is excluded 
in the summer |4*H 12 5 j . A peak is seen around 5.6 MeV/nucleon, it indicates that the 
fluctuation of the fragment distribution is the largest in this excitation energy region. 

2 

Similarly, the Normalized Variance in Z max /ZQp distribution (i.e. NVZ = ^^/Zqp) ) 
[4"2] shows a maximum in the same excitation energy region [Fig. Efd)], which illustrates 
the maximal fluctuation for the largest fragment is reached around E* / A = 5.6 MeV. 
Except the largest fragment, the second largest fragment also shows its importance in 
the above turning point. Fig. |2fe) shows a broad peak of (Z 2m ax) - the average atomic 
number of the second largest fragment exists at 5.6 MeV/nucleon. 

2.3. Fragment Hierarchical Distribution: Zipf plot and Zipf law 

In addition to the largest fluctuation, observables revealing some particular topological 
structure may also reflect the critical behavior for a finite system. In this section we 
discuss a new observable and law: nuclear Zipf plot and Zipf law. 

Zipf's law has been known as a statistical phenomenon concerning the relation 
between English words and their frequency in literature in the field of linguistics [43J. 
The law states that, when we list the words in the order of decreasing population, the 
frequency of a word is inversely proportional to its rank [33! • Recently, Ma proposed 
measurements of the fragment hierarchy distribution as a tool to search for the liquid 
gas phase transition in a finite nuclear system based on the calculations with isospin- 
dependent lattice gas model (I-LGM) and classical molecular dynamics model (231 S3]- 
The fragment hierarchy distribution can be defined by the so-called Zipf plot, i.e., a plot 
of the relationship between mean sizes of fragments which are rank-ordered in size, i.e., 
largest, second largest, etc. and their rank [2~3] . 

In such an analysis, the cluster size is employed as a variable to make a Zipf-type 
plot, and the resultant distributions can be fitted with a power law, 



M k = Akm , 



(3) 



(Z n ) oc n 



-x 



(4) 
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Figure 3. (Left part) Zipf plots for the fragments from 129 Xe at temperature T = 3 
- 7 MeV with I-LGM calculation. The histogram are simulation rseults and the lines 
are Zipf power-law fits (Eq. . From Ref. 01] . 

Figure 4. (Right part) Zipf plots in six different excitation energy bins for the QP 
formed in 40 Ar + 58 Ni. The dots are data and the lines are Zipf-law fits. The statistical 
error is smaller than the size of the symbols. From Ref. |25| . 

where n represents the rank and A is the Zipf 's law parameter. Fig. |3] shows the Zipf- 
type plots for the disassembly of 129 Xe with I-LGM calculations [44j. The fit parameter 
A was extracted in each temperature. When A ~ 1, Zipf's law is satisfied. In this case, 
the mean size of the second largest fragment is 1/2 of that of the the largest fragment; 
That of the third largest fragment is 1/3 of the largest fragment, etc. 

Now we return to the TAMU data [2^], as shown in Fig. El the region of 5-6 
MeV/nucleon excitation energy is related to the critical behavior as shown in above 
sections. The significance of this region in our data is further indicated by a Zipf's 
law analysis. Fig. H] shows the Zipf plots in six excitation energy windows. This rank 
ordering of the probability observation of fragments size of a given rank, from the largest 
to the smallest, does indeed lead to a Zipf's law power law parameter A ~ 1 in the 5-6 
MeV/nucleon range. 

In a recent analysis for multifragment emission for CERN EMU13 data, Dabrowska 
et al. analyzed emulsion data for 158 A GeV Pb-Pb and Pb-Plastic collisions and tested 
the nuclear Zipf law of Ma (2ZI- They found that their data are roughly consistent 
with the nuclear Zipf law in a certain multiplicity where the multiplicity of IMF and 
r parameter of the charge distribution reveal the same turning point. This has been 
interpretated as an evidence for the existence of the critical temperature associated with 
a liquid gas phase transition [27j . 
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Figure 5. (Left part) The deduced caloric curves for Ar quasi-projectiles (A ~ 36) in 
TAMU data. The meanings of symbols are displayed in insert. For details, see Ref. 

m- 

Figure 6. (Right part) The extraction of the critical exponent /3 (a) and 7 (b). See 
texts for details. 



2.4- Caloric Curve 

The caloric curve which relates the internal energy of an excited system at 
thermodynamic equilibrium to its temperature is a priori the simplest experimental 
tool to look for the existence of a phase transition. 

In our TAMU data the caloric curve was constructed [26J and a monotonous rising 
behavior was observed instead of a plateau-like structure observed in heavier nuclei 
The caloric curve of this lighter system looks more like cross over behavior rather than 
first order phase transition. The corresponding initial temperature is determined and it 
is 8.3 ± 0.5 MeV in the critical region of 5.6 MeV/nucleon excitation energy |26j . 

2.5. Critical Exponents 

In terms of the scaling theory, some critical exponents can also be extracted from the 
moment analysis of cluster distributions |%T) |2*5] . namely M& of Eq. El Since, for our 
system, we have already deduced a critical temperature T c = 8.3 MeV at point of 
maximal fluctuations [2E] from the caloric curve, we can use temperature as a control 
parameter to test the critical exponents. The critical exponent j3 can be extracted from 
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Table 1. Comparison of the Critical Exponents 



Exponents 3D Percolation Liquid-Gas Our Data 



T 


9 1 R 


9 91 






0.41 


0.33 


0.33±0.01 


7 


1.8 


1.23 


1.15±0.06 


a 


0.45 


0.64 


0.68±0.04 


the relation 


Zmax 


c 





(5) 



and the critical exponent 7 can be extracted from the second moment, via 

M 2 oc|l-^r. (6) 

J- c 

In each, |1 — ^-\ is the parameter which measures the distance from the critical point. 

The upper panel of Fig. |H1 explores the dependence of Z max on (1 — |r). A dramatic 
change of Z max around the critical temperature T c is observed. LGM calculations also 
predict that the slope of Z max vs T will change at the liquid gas phase transition |43j . 
Using the liquid side side points, we can deduce the critical exponent f3 by ln[Z max ) vs 
ln\l — T/T c \. An excellent fit was obtained in the region away from the critical point 
and a critical exponent (3 = 0.33 ± 0.01. Near the critical point, the finite size effects 
become stronger so that the scaling law is violated. 

To extract the critical exponent 7, we take M2 on the liquid side without Z max . 
Fig. E^b) shows ln(M 2 ) as a function of ln(| 1 — ^|). We center our fit to Eq. (0) about 
the center of the range of (1 — T/T c ) which leads to the linear fit and extraction of j3 as 
represented in Figure. EJb). The critical exponent 7 = 1.15 ± 0.06. It is seen that the 
selected region has a good power law dependence. 

Since the critical exponents (3 and 7 have been obtained, the critical exponent a can 
be deduced by the scaling relation a = -J^. In such way, o = 0.68 ± 0.04 is obtained, 
which is also very close to the expected critical exponent of a liquid gas system. 

To summarize critical exponents of our NIMROD data, we present the results in 
Table 1 as well as the values expected for the 3D percolation and liquid gas system 
universality classes. Obviously, our values for this light system with A~36 are consistent 
with the values of the liquid gas phase transition universality class rather the 3D 
percolation class. 

3. Conclusions 

In conclusion, the TAMU NIMROD data was focused and some observables relating 
to the liquid gas phase transition or critical behavior have been presented. Around 
5.6 MeV/nucleon of excitation energy for A~36 system, the fragment yield distribution 
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obeys the power law which was predicted by Fisher droplet model. From fragment size 
structure, we found that there exists a particular hierarchical arrangement, so-called 
the nuclear Zipf law, which was demonstrated by both Texas A&M heavy ion data 
and EMU13 CERN emulsion data. While, many observables demonstrate the existence 
of maximal fluctuations around this excitation energy. These fluctuation observables 
include the Campi scattering plots and the normalized variances of the distributions 
of the largest fragment (Z max ). Caloric curve was reconstructed, and then critical 
temperature and critical exponents are also discussed. Caloric curve shows a gradual 
rise of temperature versus excitation energy, no plateau is observed for such lighter 
system. Thanks to the determination of critical temperature from the caloric curve, 
critical exponents are deduced and the values are close to the liquid gas universal class. 

Finally, we should mention that a unique signal can not give any definite information 
whether the system is in a critical point or is undergoing a phase transition since we 
are facing a transient finite charged system. Only many coherent signals, such as the 
fluctuation peak, critical exponents, Fisher scaling as well as Zipf's law etc, emerging 
together can corroborate the observation of a phase transition or a critical behavior in 
finite nuclei. 
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